When is the two-level approximation untenable in issues of Decoherence? 
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We examine the conditions in favor and necessity of a realistic multileveled description of a deco- 
hering quantum system. Under these conditions approximate techniques to simplify a multileveled 
system by its first two levels is unreliable and a realistic multilevel description in the formulation of 
decoherence is unavoidable. 

In this regard, our first crucial observation is that, the validity of the two level approximation 
of a multileveled system is not controlled purely by sufficiently low temperatures. We demonstrate 
using three different environmental spectral models that the type of system-environment coupling 
and the environmental spectrum have a dominant role over the temperature. Particularly, zero 
temperature quantum fiuctuations induced by the Caldeira-Leggett type linear coordinate coupling 
can be infiuential in a wide energy range in the systems allowed transitions. The second crucial 
observation against the validity of the two level approximation is that the decoherence times being 
among the system's short time scales are found to be dominated not by the resonant but non- 
resonant processes. We demonstrate this in three stages. Firstly, our zero temperature numerical 
calculations reveal that, the calculated decoherence rates including relaxation, dephasing and leakage 
phenomena show, a linear dependence on the spectral area for all spectral models used, independent 
from the spectral shape within a large environmental spectral range compared to the quantum 
system's energies. Secondly, within the same range, the decoherence times only have a marginal 
dependence on the translations of the entire frequency spectrum. Finally, the same decoherence 
rates show strong dependence on the number of coupled levels by the system-environment coupling. 

PACS numbers: 03.65.X,85.25.C,85.25.D 



I. INTRODUCTION 



The nature of the environmental coupling and the 
properties of the environmental spectrum are essential 
factors in the determination of the decoherence prop- 
erties of multileveled quantum systems (MLS). Most 
of the realistic MLS have Caldeira-Leggett type linear 
coupling^'^ of a macroscopic system coordinate to an ex- 
ternal environmental noise field. In this publication, the 
MLS is identified by a large number of energy levels cou- 
pled by non-negligible matrix elements induced by an 
environmental perturbation. The first two levels com- 
prise the non ideal qubit space which we consider to be 
energetically well separated from the higher levels. The 
analysis of the decoherence properties of such MLS has 
recently appeared in a few publications^il. 

Normally a physical environmental spectrum (the so 
called noise) has a large range of frequencies with differ- 
ent characteristic behavior in different spectral regions. 
Moreover, the spectral noise range usually extends well 
beyond the energy scales of the quantum system that it 
interferes with. In most of the environmental models, 
weak noise signals and a realistic linear coupling suffices 
to model the system-noise interference. On the other 
hand, it is customary to consider the decohering quan- 
tum system as a two level approximated version of a 
generic MLS. The well studied spin-boson^ or central 
spini models are such examples arising from the assump- 
tion that for sufficiently low temperatures T <C Ai? a 
MLS is well approximated by its first two levels^ (where 
/S.E is a characteristic energy scale separating the first 



two levels from the higher excited states of the quan- 
tum system). This two level approximation (2LA) has 
in its essence three assumptions: a) that the incoherent 
transitions caused by the environment in the system are 
generated by the resonant processes: Implicit in the 2LA 
is the belief that the spectrum must have non-negligible 
couplings at the right transition frequencies at which the 
system makes transitions to higher levels, b) at zero 
or sufficiently low temperatures there are no available 
environmental states to couple with the system: This 
assumption and the notion of the right transition fre- 
quency basically help to neglect all parts of the spectrum 
AE < uj because of the religious belief that at sufficiently 
low temperatures the relevant part of the spectrum is the 
low energy sector w < T <C AE of which coupling is sup- 
pressed by thermal considerations, c) negligible leakage 
of the wave-function to higher levels: In fact these effects 
have not been seriously questioned until recently. Con- 
trary arguments to the negligible leakage can be found 
for instance in the recent publication Ref . [^] . 

This manuscript presents counter arguments to these 
three fundamental assumptions of the 2LA mentioned 
above by considering a MLS interacting with a bosonic 
environmental bath at zero temperature through a 
Caldeira-Leggett type linear coordinate coupling. By 
considering a spectral noise range extending well beyond 
the thermal and the system energy scales, and adopt- 
ing a model MLS for which T <C AE is well satisfied we 
demonstrate that the three basic assumptions of the 2LA 
are neither necessary nor sufficient for its validity. 

The method that we employ is the direct numeri- 
cal solution of the master equation for the multileveled 



quantum systems' reduced density matrix (RDM) in 
the short-time Born-Oppenheimer approximation. That 
the calculated decoherence rates are largely independent 
from where the energy levels are, is demonstrated. This 
observation of the independence of the decoherence rates 
from the location of the resonant transition energies pro- 
vides the first clue of the strong influence of the non- 
resonant processes. The main results are two fold: i) 
Decoherence is dominated by the non-resonant processes 
at zero temperature implying that decoherence times are 
finite and saturate near zero temperature: This observa- 
tion basically invalidates the assumption (a) above, ii) 
In a realistically multileveled system, and following (i), 
non-resonant processes induced in the entire environmen- 
tal spectral range additionally cause wave-function leak- 
age from the qubit subspace to the system's higher lev- 
els. Our numerically calculated leakage (L) times are in 
the same order of magnitude as the relaxation/dephasing 
(RD) times. Therefore the leakage, being a characteris- 
tic MLS effect, renders itself to be non-negligible at short 
times. The observation of strong leakage at short times 
invalidates the conventional belief that the affect of mul- 
tilevels in decoherence is a mere renormalizion of the ef- 
fective two level system parameters. 

The zero temperature quantum fluctuations of the 
model system in ^ are caused by the system- 
environment coupling in which the number of excitations 
in the environmental modes is unconserved by the in- 
teraction. Decoherence is caused by the participation of 
large number of these modes in the interaction. Similar 
zero temperature decoherence mechanisms have been re- 
cently verified for the mesoscopic persistent current rings 
experimentally* and studied theoretically^. In particular, 
the saturation observed in the electron dephasing time in 
disordered conductorsiS. has been argued in favour of the 
zero temperature quantum fiuctuations^^. 

In Section II. A the observation (i) is demonstrated 
starting from a finite level system- environment model. 
Most of the results are confined to three levels. Section 
II. B is devoted to the explicit multilevel effects which 
demonstrate the observation (ii). In section II. C the in- 
fluence of the non-resonant processes and the multilevels 
on the noise vacuum fluctuations arc examined. 



II. A MULTILEVEL SYSTEM-ENVIRONMENT 
MODEL 



At the first step we show the strong influence of the 
non-resonant transitions by finding that the decoherence 
times are basically independent from the shape and the 
location of the spectrum. The calculated relaxation, de- 
phasing and leakage (RDL) rates are linearly dependent 
on the spectral area in a large range of environmental 
frequencies. Therefore the participation of the system's 
higher transition levels is unavoidable. In the second 
step, we explicitly calculate the RDL rates as function 
of the number of levels in MLS and demonstrate that the 



higher transition levels can strongly affect the rates. 

Although some results are presented for M leveled 
MLS, most of our conclusions are based on the solution 
of the three level system-environment Hamiltonian (i.e. 
M = 3) 
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Here Ei — E2 ~ 0.5 and E^ — 2.5 are fixed system 
eigenenergies in some absolute (and irrelevant) energy 
scale and |r) indicates the r'th system eigenstate. We 
consider the couplings fnr to be real-symmetric and 
¥^12 = ^23 = 0.1, and (pi3 = in the same absolute en- 
ergy unit. The frequency dependent part of the system- 
environment coupling is separately denoted by rji^. All 
other energy (time) scales -including the spectral en- 
ergy scales are dimensionless in the same absolute en- 
ergy (time) scale. The environment is characterized by 
a spectral function I{uj) — rj^ (2n(^ -I- 1) where n^^ is the 
bosonic thermal distribution. Since our results are con- 
fined to zero temperature the bare distribution n^ — Q, 
and therefore I(oj) = rj^ for which we use two spectral 
models below. For obtaining results in regard of the pure 
two level system we simply reset the couplings to the 
third level in ^ to zero (i.e. (^13 = (^23 — 0). From the 
2LA point of view, we have a quantum system (Q which 
is supposed to be a good candidate for a two level system 
(Since AS — E3 — E2 j^ and the environmental tem- 
perature is zero; i.e. T ^ AE is well satisfied). However, 
as we show below, the model ^ breaks all three aspects 
of the 2LA. 

The machinery that we use is the master equation for- 
malism in the interaction picture. The quantum system 
is assumed to be initially prepared in the qubit subspace 
IV'(O)) = a|l) -I- 6|2) where we consider a = -s/OT , 6 ~ 
e*'^' ^ VO.O. For the short time dynamics it is sufficient to 

use the Born-Oppenheimer approximationpiS in the time 

- (S) 
evolution of the density matrix p{t) — p (t) (g) /5g(0). 

Here " denotes the interaction picture, p^"^^ denotes the 

system's reduced density matrix (RDM) and pe is the 

density matrix of the pure environment characterized by 

the spectral function /(oj). The master equation for the 

system's RDM is^ 






dt' 



r,s 



Kr(t,t')piV{t') (2) 



where 

= \T(t - t')[((ft'ft')nr^s,m ^ ('ft')nr('ft)sm] (3) 

+J^*{t - t')[{(ft'<ft)msSr,n " i'ft)nri'ft')srn]j 

is the non-Markovian system-noise kernel. Here {<pt)nr — 
(<^o)nr e"*'^""'^''-'*. The noise correlator J-{t) = 



0(t) J_^ du}e^'^^I{uj) is calculated for three different 
types of spectral functions: a) a rectangular spectrum 



and b) a Lorentzian spectrum 



{UJ - LUo)^) 



.{uj) = 



A 



n {oj - UJq)'^ + £2 



and the power-Gaussian spectrum 



m^Ai^^r-^ 



,-u;V4A2 



A' 



(4) 



(5) 



(6) 



The common parameters in the first two spectra are the 
location of the center of the spectra coq and the spectral 
width e. The spectral area in both cases is denoted by 
A = J_ dujlluj). The power-Gaussian spectrum is sim- 
ilar to the original model spectrum studied by Caldeira 
and Leggett. This spectrum, along with the power ex- 
ponential, is one of the most frequently used spectral 
models nowadays. The spectral area in ^ can be cal- 
culated analytically for different values of the Ohmicity 
parameter v and the Gaussian cut-off A. The normaliza- 
tion with respect to the spectral area is included in the 
overall constant A. 

The relaxation, dephasing and the leakage rates are 
calculated by examining P22 (t),\p\2 {t)\ and 1 — [p^ ^ ^ (i) -I- 
P22 {t)] respectively at short times 1 < t < 2.5 [in units of 
the inverse absolute energy scale in (^] where the domi- 
nant behavior is well-fitted to exponential. 



A. The influence of the spectrum on RDL 

The decoherence rates (t^^^)"^ corresponding to the 
qubit subspace of the system found by the numerical so- 
lution of ||2J) are studied below. In Fig.J^I the rates are 
given for a varying spectral center cjq with fixed spectral 
areas. The figure suggests that the decoherence rates 
in the qubit subspace are largely independent from the 
center of the spectrum for the model spectra in Eq's.lQ)) 
and © and for the indicated range of spectral areas. 
In Fig.® and for the same three level system the RDL 
rates are plotted against the spectral area for the rectan- 
gular spectrum. The curves obtained for the Lorentzian 
spectrum totally coincide with the corresponding ones for 
the rectangular spectrum and they are not shown. On 
the other hand, the same quantities are also plotted in 
Fig. (PJ for the data obtained from the power-Gaussian 
spectrum for u — —1,0,1. The shown relaxation rates in 
Fig. © are indistinguishable from those in Fig. (|2J) within 
a very large range of spectral areas (The calculated de- 
phasing and the leakage rates are in the same order of 
magnitude as the relaxation rates which can also be fol- 
lowed from Fig. (|2l; hence, they are not shown in the 
figure) . These observations prove the independence from 
the shape of the spectrum. Moreover, the plots (a) for 
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FIG. 1: Dependence of the RDL rates on the spectral shifts 
for various spectral areas. Note that the vertical axis is loga- 
rithmic. 



ly — 1, (h) for z^ = and (c) for i/ = — 1 in Fig.Q yield 
identical decoherence rates for the same horizontal scale. 
Hence, in the studied range — 1 < i' < 1 corresponding 
to sub-Ohmic, Ohmic and super-Ohmic (i.e. j^ = — 1, 0, 1 
respectively) the decoherence rates are determined solely 
by the spectral area. This observation is contrary to the 
earlier belief that it is the low frequency sector in the 
environmental spectrum [dominated by the w^"*"" term in 
©] dominating the RDL rates. 

We hence observe independence from the shape and 
the type of spectrum within a large range of spectral ar- 
eas covering three orders of magnitude. The log-log axes 
in Fig's (O and ^ additionally indicate that the RDL 
rates depend linearly on the spectral area. The mono- 
tonic dependence on the spectral area (whether it is a 
linear dependence or not), is a strong signature for the 
dominance of the non-resonant processes in the decoher- 
ence rates. Additionally, in Fig.®, and also implicitly 
in Fig.Q, the leakage rates are shown to be in the same 
order of magnitude as the RD rates. The leakage rates 
increase by the inclusion of higher levels [also see Fig. Q) 
in Section II. B below]. Hence, we conclude that, the leak- 
age demonstrates itself to be a non-negligible short time 
effect manifested in multileveled systems. 

The dominance of the non-resonant processes and their 
survival at zero temperature implies that the influence 
of the higher levels cannot be avoided independent from 
how well the qubit is energetically separated from those 
levels. This observation is in contrast with the assump- 
tion that, in low temperatures T compared to the sepa- 
ration AE between the qubit subspace and higher states 
(i.e. T ^ AE), the two level approximation is well sat- 
isfied. 
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FIG. 2: Dependence of the RDL rates on the spectral area for 
the spectral function @. Note that the both axes are loga- 
rithmic. In generating the ten level results we assumed that 
the system's energies above the first two levels are harmonic 
and given by En = (n — 1/2) , 3 < n and all the couplings are 
(finr =0.1 for n + r = odd and zero otherwise. 
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FIG. 3: Dependence of the relaxation rates on the spectral 
area for the spectral function ^ corresponding to a) ;/ = 1 
super-Ohmic, b) :/ = Ohmic, c) ly — —1 sub-Ohmic spectra 
and for two, three and ten level systems. The trend for the 
dephasing and the leakage rates are very similar to Fig.ll^J 
and they are not shown here. Note that the both axes are 
logarithmic. The energies E„ and couplings ipn r of the multi- 
level system are as indicated in Fig.ll^J. The horizontal range 
in (c) is common to all graphs. 



FIG. 4: RDL rates at T = against the number of levels 
for different spectral areas and different multilevel coupling 
ranges A. Here the coupling matrix elements between the 
levels n and r are modelled by (p„r = 0.1 e~'"~'~'' for 
n + r = odd an zero otherwise. Note the logarithmic verti- 
cal axis. The open and solid symbols superimposed on the 
solid line refer to dephasing and relaxation rates respectively. 
The pure symbols without connecting lines refer to leakage 
rates. Also {A, A) values are (5, 8) for squares, (5, 40) for dia- 
monds, (50, 8) for circles and (50, 40) for triangles. Here, the 
Lorentzian spectrum is used. 



B. The influence of the multilevels on RDL 

The influence of the multilevels is examined by using 
a model coupling matrix (^ in Q as 



^7, 



^O.le^l"-''!/'^ 



(7) 



for n + r — odd and zero otherwise. By varying the multi- 
level coupling range A, the effect of the actively coupled 
levels on the RDL rates can be examined. In Fig.Q 
the data is generated for different values of the multi- 
level range A — 8,40 and the spectral area A = 5,50. 
The first observation is that, the linear trend observed 
in Fig. 101 against the spectral area is respected within 
a larger range of multilevels. Secondly, the multilevel 
coupling range is observed to give rise to a saturation in 
the RDL rates at the onset M ~ A. Fig.Q) therefore 
demonstrates further evidence of the observable effects 
in the decoherence rates arising from the finite coupling 
to multilevels. 



C. The fluctuations in the noise vacuum 

A characteristic feature of the coordinate coupling is to 
produce weak vacuum fluctuations as well as finite num- 
ber of environmental modes at zero temperature. These 
corrections, can be calculated perturbatively and they 



give a clue about the effectiveness of tfie virtual (non- 
resonant) processes. 

The corrections to environmental number of modes can 
be calculated using the retarded Greens function 

GUt - t') = -iQ{t - t') {Tal{t)hUt')S{(x,, -(x))) (8) 

where T denotes the time ordering, the "" denotes the 
interaction picture and the S matrix includes the inter- 
action Hamiltonian in Q in the interaction picture. The 
environmental photon number is than found by the stan- 
dard method n^ = ilunt'^t Guj(t — t'). The average is 
taken over the noninteracting initial state including sys- 
tem and the environment in the product form 



a) 



fall) 



■b\2))® \env) 



(9) 



where \env) is characterized by the spectral function I{ijj) 
at zero temperature. Evaluating the time ordered inte- 
grals arising from the expansion of the S matrix, and 
taking the limit t' -^ t we obtain a second order correc- 
tion as 



.(2) _ 



21{ui 
A? 



X i\a\'^(pos + |6p(/?is + {a*b + b*a)ipis(pos) 



(10) 



where A^ ^ uj + Ei — Es and and Ei — E2 as announced 
in (P). Eg. ltTUI) yields a divergent result for the total 
number of photons even for arbitrarily small system- 
environment couplings. The divergence is obtained inde- 
pendently from the type of the spectrum in Q or |(SJ and 
arises from the second order singularity in H10|l . To avoid 
this unphysical result in the second order, we convert the 
second order correction into an RPA sum following con- 
ventional practices in perturbative approaches'^. As we 
use RPA, the second order pole in (|10|l is split into two 
first order poles which yields a finite result at the RPA 
level. The relevant RPA graphs are depicted in Fig. (O 
below. 

Next, we calculate the fluctuations in the number of 
photons An„ by using a four point retarded Green's func- 
tion for the environment as 

Du^it^t') = -iQ{t-t') {Tal{t)dUt)al{t')ak{t')S{^,-^)) 

(11) 
The connected part of (|ll|l directly yields the fluctuations 

in the photon number as AnJ = i limt/^t D^{t — t'). The 
diagrams corresponding to the RPA scheme using 111|) 
are shown in Fig.©. The fluctuation in the total photon 
number is found by (AA^)^ — J dtu AnJ and is given by 



(ANf 



E 



duj 



'-00 [^ + As) 1 - 4:{system)J{Lj)/{u; + As) 

(12) 
which is plotted in Fig.@ with respect to e and param- 
eterized by loq for the two, three and ten level systems. 
The bare spectral area is fixed at A = Stt. The observed 
independence from e at fixed spectral area implies inde- 
pendence from the spectral shape. Additionally indepen- 
dence from ujQ is also shown. 



where 
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FIG. 5: Contributing diagrams in the calculation of Ijll^ . (a) 
The definition of the basic second order root diagrams that 
are used to implement RPA. The dashed arc lines indicate the 
retarded limit t' ^ t + 0~ in the calculation of ni, , (b) the 
RPA hierarchy using the root diagram in (a) , (c) the resulting 
RPA Green's function used in the calculation of 1121 . Here the 
dashed arc lines indicate again the retarded limit t' ^ t + 0~ 
in the calculation of nL 
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FIG. 6: Fluctuations in the total photon number versus the 
spectral width for various ujo values computed for two three 
and ten leveled MLS. Here the bare spectral area A = 5n 



In Fig.Q below, the ratio AN/N is plotted as a func- 
tion of the spectral area. The figure, together with the 
inlet for A7V versus spectral area, is well fitted to 



AN ^ Const. 



A 



(13) 



We observe that AA^ ex V^ and A^ ex A in a large range 
of spectral areas. Here the positive constant Const, de- 
pends on the number of multi levels. 

The results of this subsection are in a sense also a check 




levels. The presented results are obtained from the nu- 
merical solution of the master equation with the Born- 
Oppenheimer approximation, and hence, they are exact 
in the short time limit where the dccohcrence properties 
are studied. 
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FIG. 7; The ratio AN/N versus the spectral area (inlet: AiV 
versus the spectral area). 



for the validity of the Born-Oppenheimer approximation. 
The correction obtained to the spectral area as a result of 
the system-environment coupling is bounded by a few ten 
percent. However, as the number of levels increase the 
total fluctuations decrease relative to the increase in the 
total number of environmental modes. For a MLS with 
strong couplings between a large number of levels the va- 
lidity of the Born-Oppenheimer approximation may still 
survive although we did not check this result explicitly. 



III. CONCLUSIONS 



It is shown that for the most conventionally used 
system-environment couplings of linear coordinate type 
the non-resonant processes overwhelm the resonant ones 
in their contribution to the decoherence rates at zero tem- 
perature. In this regard, and within the studied spectral 
range, it is observed that the decoherence rates do not 
depend on the specific low/high energy properties of the 
spectrum, whereas, a strong (and dominantly) depen- 
dence is observed on the overall spectral area. These 
results are confirmed for three independent spectral pro- 
files. 

We have also examined the effect of the system- 
environment coupling. We observed that the number 
of levels which are coupled by the system-environment 
coupling plays a non-negligible role in decoherence time 
rates. We see that as the number of coupled levels in 
the system increase, both the decoherence time rates 
and photon number corrections increase further from the 
rates and corrections in a pure two level system. The 
observed effects for the MLS cannot be explained by an 
equivalent renormalization of the two levelled subsystem. 
In the light of these observations, we conclude that the 
conventional postulates of the two-level approximation 
do not necessarily lead to a dynamical behaviour of a 
multilevelled system largely confined to its lowest two 
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